Strong Griffiths singularities in random systems and their relation to extreme value 

statistics 
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We consider interacting many particle systems with quenched disorder having strong Griffiths 
singularities, which are characterized by the dynamical exponent, z, such as random quantum 
systems and exclusion processes. In several d — 1 and d — 2 dimensional problems we have calculated 
the inverse time-scales, r , in finite samples of linear size, L, either exactly or numerically. In all 
cases, having a discrete symmetry, the distribution function, P(t ,L), is found to depend on the 
variable, u = r~ 1 L Z//d , and to be universal given by the limit distribution of extremes of independent 
and identically distributed random numbers. This finding is explained in the framework of a strong 
disorder renormalization group approach when, after fast degrees of freedom are decimated out the 
system is transformed into a set of non-interacting localized excitations. The Frechet distribution of 
P(r~ , L) is expected to hold for all random systems having a strong disorder fixed point, in which 
the Griffiths singularities are dominated by disorder fluctuations. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 

In interacting many-particle systems quenched (i.e. 
time independent) disorder can induce unusual physical 
properties which do not exist in non-random systems. 
For example in a random ferromagnet the linear suscep- 
tibility, x, is divergent in an extended part of the para- 
magnetic phase, which is called Griffiths phaseQ. This 
type of singular behavior is caused by rare regions of the 
sample, in which due to strong disorder fluctuations the 
local couplings are much stronger than their average and 
they are even stronger than the disordering field, such 
as the temperature in a classical system. Consequently 
these rare regions are locally in the ferromagnetic phase 
and the relaxation process associated to them involves 
a very large relaxation time, r. In the thermodynamic 
limit the distribution of the r-s has no upper limit which 
leads to singularities in several average physical quanti- 
ties. One can say that these Griffiths singularities are 
controlled by a line of semi-critical fixed points, having a 
diverging relaxation time but a finite correlation length, 
£. In a classical disordered system in which the phase 
transition is triggered by the variation of the tempera- 
ture these singularities are very weakQ and often one 
needs astronomical times to be able to observe them. 
There are, however, another type of random systems in 
which the Griffiths effects are much stronger and one can 
measure them in conventional experiments. 

Strong Griffiths effects can be observed, among oth- 
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ers, in random quantum systems^ at zero temperature, 
T = 0, in which a quantum phase transition can take 
place by varying a control parameter, S. The prototype of 
random quantum systems is the random transverse-field 
Ising model (RTIM), which is defined in one dimension 
(Id) by the Hamiltonian, 
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in terms of the Pauli matrices, cr^ ,z , at site i. The ex- 
change couplings, Aj, and the transverse fields, hi, are 
both independent and identically distributed (ii(i)ran- 
dom variables and the control parameter is defined 4] as: 



5 = 



[m/t]av - [lnA] av 
var[ln A] + var[ln h] 



(2) 



Here and in the following we use [. . . ] av to denote aver- 
aging over quenched disorder and var[:c] stands for the 
variance of x. The system is in the ferromagnetic (para- 
magnetic) phase for <5 < (<5 > 0) and the quantum 
critical point is located at 8 = 0. It was first McCoy[5| 
who calculated exactly the average susceptibility of this 
model and found that [x]av is divergent in a finite regime 
of the paramagnetic phase. Later the excitation energy, 
e ~ f/r, is shown 4] to vanish with the size of the system 
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with the dynamical exponent, z > 0, which is a continu- 
ous function of 6. 

Another class of problems with strong Griffiths effects 
are stochastic many-particle systems Q) with quenched 
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disorder 7], such as the Id partially asymmetric sim- 
ple exclusion process (PASEP)) with position or par- 
ticle dependent hopping ra tes [a or the zero range 
process (ZRP) with disorder[l(j. Also reaction-diffusion 
models|ll| with quenched disorder mig ht show strong 
Griffiths effects 121, see, however [TsL 141. As an exam- 
ple we consider here the PASEP with particle-wise (j>w) 
disorder 8], in which the i-the particle, % = 1, 2, . . . , N, 
can hop to empty neighboring sites with a rate pi to the 
right and to the left in a Id periodic lattice with L > N 
sites. Here the hop rates are iid random numbers and the 
control parameter is defined as: 

s _ [mp] av - [lnq] av 
p var[lnp] + var[ln q] 

so that the particles move to the right (to the left) for 
S p > (Sp < 0) . Here we restrict ourselves to the domain, 
S p > 0. Strong Griffiths effects in the PASEP are caused 
by such rare particle clusters in which the local hopping 
rates of the particles prefer jumps against the global bias. 
As a result in the anomalous diffusion regime, which can 
be called as the Griffiths phase of the PASEP, the average 
displacement of particles grows in time as 

(x(t)) ^t 1/{1+z "\ (5) 

when the system undergoes a coarsening process. Here, 
the dynamical exponent of the PASEP, z p > 0, is a con- 
tinuous function of S p . Note that outside the Griffiths 
phase, S > Sg, the stationary velocity of the particles is 
v = 0(1), whereas in the Griffiths phase it vanishes asQ: 

v ~ L" z » . (6) 

In general Griffiths singularities are characterized by 
the distribution of the relaxation times of slow processes 
or equivalently by the distribution of their inverse, 1/r, 
which is just the excitation energy, e, of the RTIM and 
the stationary velocity, v, of the PASEP. In a finite sys- 
tem 1/t is just the smallest event associated to one of 
the rare regions and the occurrence of these rare regions 
as well as the distribution of 1 /r depends on the form of 
the randomness in the dynamical model. In the mathe- 
matical literature such questions are studied in the frame 
of extreme value statistics [l^] (EVS) and rigorous results 
are known for iid random variables, which are distributed 
according to a parent distribution, n(y) (y = 1/t). In 
this case the limit distributions are in a few different 
universal scaling forms depending on the asymptotic be- 
havior of the parent distribution for large values of y. If 
the decay of 7r(y) is faster than any power then it is the 
Gumbel distribution, for power-law tails it is the Frechet 
distribution and for power-laws with an edge the Weibull 
distribution is the result. If the random variables are 
not iid there are no general mathematical results, how- 
ever a few special cases have been studied and solved [l5j. 
Among others we can mention the statistics of extreme 
intensities of Gaussian interfaces [ljj which generally does 



not coincide with the known limit distributions of iid vari- 
ables. Similarly, for hierarchically correlated variables 
the limit distribution deviates from the Gumbel formal • 

As we described above strong Griffiths affects are 
present in interacting particle systems having strong tem- 
poral and spatial correlations, therefore at first thought 
it seems unlikely that a simple connection with the the- 
ory of EVS exists. A closer look, however, shows that in 
a class of systems the rare events, which are the source 
of Griffiths singularities, are well localized and separated 
in such a way that they could be considered indepen- 
dent. For these systems then the distribution function, 
Pl (1/t), could correspond to the appropriate limit dis- 
tribution of EVS. 

In this paper we are going to study this issue in more 
details. For some random Id models having a small corre- 
lation length we perform exact calculations. On the other 
hand at a general point of the Griffiths phase a renor- 
malization group (RG) approach [18j is applied, which is 
expected to give asymptotically exact dynamical singu- 
larities in the Griffiths phase. The RG equations can be 
analytically solved for some Id random models, such as 
for the RFIM0 and for the PASEP @, in other cases 
and in higher dimensions one resort to numerical calcu- 
lations. For some Id random models we also calculate 
numerically the different distribution functions. 

The structure of the paper is the following. In SeclTTl 
we consider Griffiths singularities in random quantum 
systems. In more details we study the Id RTIM, but also 
the random quantum Potts chain, as well as the RTIM for 
ladders and in 2d and the random Heisenberg model in 
Id and 2d are investigated. In Sec lllll we study Griffiths 
singularities in random stochastic systems, such as in the 
PASEP with particle-wise and site-wise disorder. Our 
paper is closed by a discussion in Sec II VI 



II. GRIFFITHS SINGULARITIES IN RANDOM 
QUANTUM SYSTEMS 

A. The RTIM in Id 

1. Exact results for extreme disorder 

In order to gain some insight into the origin of Griffiths 
singularities we solve exactly the RTIM in Eq.QJ using 
the bimodal distribution: Aj = A, with probability c and 
Aj = A -1 with probability 1 — c, and hi = 1. Furthermore 
we take the limit c«l, when most of the couplings are 
weak and A 3> 1, so that we are deeply in the paramag- 
netic phase. There are, however, rare regions in which 
all the couplings, say n, are strong and appear with an 
exponentially small density, p(n) — c™. Such a cluster 
is typically embedded into a see of very weak couplings, 
thus the corresponding excitation energy is obtained by 
solving the gap in a cluster with free boundary condi- 
tions, leading to: e(n) w A - ™. Now we obtain that the 
distribution of the excitation energies in an infinite sys- 
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tem has a power-law tail: 

and the gap exponent, u>, is given by: 
_ ln(l/c) 



In A 



- 1 



(7) 



(8) 



In a finite system of size, L, the typical size of the largest 
cluster, 7ii, is given by the condition, L^2 n>n P( n ) = ^> 
so that the typical value of the smallest gap, ei, is given 
by: t\ = L~ z , with the dynamical exponent defined in 
Eq.@: 



In A 



(9) 



ln(l/c) 

Note that we obtain the relation, uj = 1/z — 1, so that 
the distribution function for finite systems satisfies the 
scaling form[l^ : 



P L (e 1 ) = L 2 P 1 (eii 2 ), 



(10) 



and the subscript of the scaling function, Pi(u), refers 
to the first gap. In the following we determine the scal- 
ing function. First, we note that the rare regions in the 
system are localized thus they could be placed at ~ L dif- 
ferent positions of the chain. Furthermore the different 
rare regions are independent and their lengths, which are 
proportional to the log-gap, is distributed by the same 
exponential distribution. From these follows that t\ is 
the smallest gap out of ~ L independent possible gaps, 
associated to the different rare regions, and the gaps are 
distributed by the same parent distribution given in a 
power-law form in Eq.Q. Consequently in the large- 
ly limit the distribution function, Pi (it), is given by the 
standard Frechet distribution[15|: 



Pi(u) = i«V*-i eX p(-u 

z 



1M 



(11) 



and u = uq€iL z , where the non-universal constant, uq, 
depends on the amplitude of the tail in Eq.JTJ). 



2. Analytical results of the strong disorder RG 

The strong disorder RG method^ has been intro- 
duced by Ma and Dasgupta^ to study random an- 
tiferromagnetic Heisenberg chains. For the Id RTIM 
FisherQ has applied the method and solved analytically 
the RG equations in the vicinity of the random quantum 
critical point. During renormalization the strongest local 
term in the Hamiltonian, coupling or transverse field, is 
decimated out and between remaining sites new interac- 
tions are generated through perturbation calculation. If 
a strong bond, say A2, connecting sites 1 and 2 is deci- 
mated out the effective transverse field, h, acting on the 
two-site cluster is given by 4]: 



hih 2 

K\2 



(12) 



where k — 1 . On the contrary if a strong transverse field, 
say hi and thus the site 2 is decimated out the effective 
coupling, A, connecting the remaining sites, 1 and 3 is 
given by 3: 



A = 



A1A2 
nh 2 



(13) 



with k — 1 . As renormalization goes on the energy scale, 
set by the strength of the strongest local term in the 
Hamiltonian, gradually decreases and the distributions 
of the renormalized couplings and transverse fields ap- 
proach their fixed-point form. At the critical point cou- 
plings and transverse fields are decimated symmetrically 
and their fixed-points distribution functions are identical 
and have the so called infinite disorder property 3 : the 
ratio of any two terms tends either to zero or to infinity. 
Consequently the decimation rules become exact and the 
critical singularities, both static and dynamical ones are 
asymptotically exact. In the disordered Griffiths phase, 
which is a part of the paramagnetic phase after a start- 
ing period of the RG when spin clusters of typical size, £, 
are formed almost exclusively transverse fields are deci- 
mated out. At the fixed point the energy scale still goes 
to zero and the decimation equations are asymptotically 
exact leading to exact dynamical singularities, however 
the spatial correlations are short-ranged and correct only 
up to a range of £. The analytical solution of the RG 
equations by FisherQ is exact up to 0(S) which is then 
extended into the complete Griffiths phase |2lL |22|. The 
dynamical exponent, 2, is found to be the positive root 
of the equation j2^: 



l/z 



= 1 



(14) 



which has also been derived by a mapping with the ran- 
dom random walk|23j . For the binary disorder used in the 
previous section this leads to: A 1 ^ z c + A _1 / Z (l — c) = 1, 
having the solution in Eq.@ in the limit c <C 1. Note 
that this result is valid for any value of A > 1, which 
is connected to the fact that the gap of a large n-bond 
cluster in a see of weaker bonds scales as: e(n) ~ A - ™. 
The dynamical exponent in Eq. (|14|l is a continuous func- 
tion of 5, and in the vicinity of the critical point it is 
divergent as 4] z m 1/25. It can be shown by the RG 
calculation that the low temperature behavior of the av- 
erage susceptibility and that of the specific heat, cy, are 
given bvplE^: 



X(T) ~T- 1+1 / 2 , cy^-T 1 / 2 



(15) 



thus the singularity involves the dynamical exponent. 
Note that at zero temperature x(0) is divergent for z > 1. 

The distribution of the smallest energy gap in a finite 
chain of length, L>^, and in the vicinity of the critical 
point: < S ^ 1 is calculated by Fisher and Young |24|. 
For the distribution of the log-gap, G = — lnei, they 
obtained: 



dProb[G|L] w 2SLn G e^q>(~Ln G )dG , 



(16) 
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where 

n G ks 4(5 2 exp(-2<5G) (17) 

is the density of remaining clusters. Here making use 
of the fact that in the vicinity of the critical point the 
correlation length is given by[J £ ~ 8~ 2 and z w 1/25 
we obtain: 



Inserting Eq. l|18|l into Ea. i|16|) we obtain that the scal- 
ing variable is u = uq6i(L/^) z , in terms of which the 
distribution is in the Frechet form in Ea. i|ll|) . 

3. Phenomenological considerations 

Results of the previous two subsections indicate that 
the distribution of the smallest gap of the RTIM in Id 
is the Frechet distribution, at least in the limiting cases 
i) extreme binary disorder with £ <C 1 but z is arbi- 
trary fSec lII A~T|) ii) arbitrary form of disorder, but 6 <C 1 
fSec lII A2~\i . Here we argue that the Frechet form of the 
distribution should be generally valid in the disordered 
Griffiths phase, at least if the low-energy excitations in 
the system are localized. Here we note that performing 
the RG up to an energy scale, Oo, which is much smaller 
than ~ £~ z , we obtain an equivalent random spin 
chain with very weakly interacting effective spins. Typ- 

—Xlz 

ically, the log-bonds are in the order of — In A ~ fl Q , 
so t hey are small compared to the value of effective fields 
h^lo [23] ■ Consequently, in the following RG steps almost 
exclusively transverse-fields are decimated out through 
the transformation in Eci. (|13|) . which does not influence 
the transverse fields of the other active spins. At the en- 
ergy scale, f2o> the effective transverse fields, h, have a 
power-law distribution P(h) ~ h~ 1+1 ' z , and the small- 
est one gives thus the smallest gap of the chain. Con- 
sequently the renormalized chain satisfies the conditions 
needed for the validity of the Frechet limit distribution. 
Generally, for the n-th excitation energy the limit distri- 
bution is given bv|l5j: 

P„K) = exp(-uV*) , (19) 

in terms of u n = uqL/ z c n . We note that the gap exponent 
is now, u) n = n/z—1, which has been obtained before [27] 
through scaling considerations. 

This type of reasoning for the smallest gap applies to 
another random systems, too, provided the decimation 
rules are analogous to those in Eas.l|12|l and (|13fl . In d- 
dimensions the only difference is that z is replaced by 
z/d, provided the strong disorder RG approach leads to 
localized, non-interacting effective degrees of freedom. 

Next, the distribution of excitation energies is analyzed 
in the ordered phase of the RTIM. In Id due to duality 



the disordered and ordered phases are related to each 
other and in the RG procedure the role of the couplings 
and the transverse-fields are exchanged. Consequently, 
at the energy-scale, <S Q^, the typical log-fields are 
— In h ~ fl Q , whereas A ~ S2. Thus the renormal- 
ized chain is a classical random bond Ising chain with 
very small effective transverse fields having an exponen- 
tially vanishing first gap. The dynamical properties of 
the chain are governed by kink- like excitations [25ij. In an 
open chain the lowest excitation is one kink, which cor- 
responds through duality to a spin flip in the disordered 
phase. Consequently the distribution of the second gap 
in an open chain follows the Frechet distribution. On the 
contrary for a closed chain the lowest excitation involves 
two kinks and the corresponding excitation energy is the 
sum of the first two smallest effective couplings of the 
renormalized chain. Consequently the distribution of the 
second gap in a closed chain is related to the solution 
of the above extreme value problem. In higher dimen- 
sions the singularities in the ordered Griffiths phase have 
a more complicated structure 28] . Here the renormalized 
system is a random bond Ising model with a complicated 
topology, and the (second) excitation, which is relevant in 
the dynamical properties corresponds to the creation of 
an oppositely magnetized domain. The excitation energy 
is just the sum of the effective couplings at the boundary 
of the domain. 

In the following we check the above conjectures by nu- 
merical calculations. 



4- Distribution of gaps from numerical diagonalization 

The RTIM in Id can be transformed into a problem 
of free fermions and the calculation of the gaps necessi- 
ties the diagonalization of a 2L x 2L tridiagonal matrix 
the entries of which are the couplings and the transverse 
fields pt|. In the numerical calculation we used a contin- 
uous, uniform distribution: tt\(\) = 1 for < A < 1 and 
otherwise; as well as TVh(h) = 1/ho, for < h < ho and 
otherwise. The critical point is located at ho = 1, the 
disordered Griffiths-phase is in the region 1 < ho < 00, 
and the dynamical exponent is given by the solution of 
the equation: zln(l — z~ 2 ) = — \nho. We have calculated 
the first gaps for finite systems with L = 64, 128 and 256 
and for two values of ho — 2 and 3. The probability dis- 
tribution of the first gaps are shown in Fig^which all fit 
well to the Frechet distribution. 

The same conclusion is obtained with the distribution 
of the second gap, £2, which is presented in Fig 

B. Numerical RG study of chains, ladders and 
higher dimensional systems 

For more general models, in particular for a topology 
which is more complex than a linear chain to calculate 
the gaps one resorts to numerical implementation of the 



5 



x 



- h =2 




i 

- 




















'h =3 




h 






















■ 



-6 -4-2 2 



x=-ln(e 1 L i 



FIG. 1: Scaling plot of the distribution of the first gap of the 
Id RTIM for different sizes, L = 64 (+,□), 128 (x,0) and 
256 (*, ■) at two values of the uniform disorder. The full lines 
are the Frechet distributions having the exact values of the 
dynamical exponent, z. 
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FIG. 2: The same as in Fig^for the second gap. 



strong disorder RG method. This type of method is first 
used in Ref.|28| and implemented for a finite system in 
Ref. (2^1 ■ In the finite lattice version we use here the dec- 
imation procedure is performed up to the last spin and 
the gap of the system is identified by the effective trans- 
verse field acting on that spin. By this procedure the 
first few RG steps are approximative, which will influence 
the value of the non-universal constant, wo, but the later 
transformation steps close to the fixed point are presum- 
ably asymptotically exact. In the numerical calculations 
we used the uniform distribution and considered at least 
10000 independent realizations of disorder. 



1. Random quantum Potts chain 

A simple generalization of the RTIM for g-state spin 
variables: \si) = |1) = |2) = . . . \q) is the random quan- 
tum Potts model defined in Id by the Hamiltonian: 



L-l 



l , q- 

hi 



1 q k - 



(20) 



where: Mi\si) = |sj + l,mod q). The control parame- 
ter of the model is in the same form as for the RTIM in 
Eq. |2j| . The strong disorder RG approach is used for this 
model in Ref. [3(| for the quantum critical point whereas 
in Ref.[2l| for the Griffiths phase. The transformation 
rules for the bonds and external fields are of the form 
given in Eqs.^l an( I l|12(l with k = 2/q, thus the phe- 
nomenological argumentation directly apply here. The 
only difference comparing with the RTIM is the degener- 
acy of the excited states. The distribution of the first gap 
is thus still of Frechet type, and gaps between subsequent 
multiplets behave as the higher excitations of the RTIM. 
For the distribution of the gaps in the ordered Griffiths 
phase one can obtain similar conclusions as described in 

sec irroi 

In the numerical application of the RG technique we 
have calculated the distribution of the gaps at a large 
finite system, L = 2048, for the uniform distribution 
with h = 3 but for different values |3lT] of q = 2,3,4 
and 8. The dynamical exponent is q dependent and 
calculated by a numerical integration of the analytical 
RG equations [32|. In Fig|3]the numerically obtained gap 
distributions are compared with the Frechet distribution 
and excellent agreement is found. 



ST 




FIG. 3: Distribution of the log-gaps, identified as the trans- 
verse field of the last remaining spin in the RG procedure, for 
the the random quantum Potts chain with different values of 
q having a uniform distribution with ho — 3.. The Frechet 
distributions are indicated by full lines for which the dynami- 
cal exponent is obtained by the solution of the analytical RG 
equations. 



6 



2. RTIM: ladder and 2d system 

Here we consider the RTIM with a more complex topol- 
ogy, first a ladder composed of two chains and afterwards 
a large plaquette of a square lattice. 

In FigQ] the distribution functions of the log-gaps of 
the ladder model with h a = 2.5 are presented in a log-log 
scale for different lengths up to L = 1024. The curves for 
different L-s are shifted to each other and a good scaling 
collapse can be obtained by using the scaling combina- 
tion in Ea. (|l()fl with a dynamical exponent, z — 2.9, as 
illustrated in the inset of FigQ] The (absolute value of 
the) asymptotic slope of the curves for small gaps is given 
by io + 1 , which for localized excitations should be related 
to the dynamical exponent as: 

d 

w + l = -, (21) 
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FIG. 5: (Color online) The same as in Fig^] for the square 
lattice RTIM with ho = 9. and for different sizes. In the inset 
the optimal scaling collapse is obtained with z = 2.7. The 
asymptotic slope of the curves is given by u + 1 = 0.74, so 
that the relation in Eq, l2H is well satisfied. 



for a d-dimensional system. Here we obtained u) + 1 = 
0.34 so that the relation in Ea. l|21() with d = 1 is very 
well satisfied. Finally, the scaled curve in the inset of 
Figgis very well described by the Frechet distribution. 
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FIG. 4: (Color online) Distribution of the log-gaps, identi- 
fied as the transverse field of the last remaining spin in the 
RG procedure, for the ladder model with ho = 2.5 and for 
different lengths. In the inset a scaling collapse according to 
Eg. 1101 with z = 2.9 is presented. The Frechet distributions 
are indicated by full lines. The asymptotic slope of the curves 
is indicated by a dashed line and is given by u + 1 = 0.34, so 
that the relation in Eq, l2H is well satisfied. 

Results of the same analysis of the gaps of the square- 
lattice RTIM for ho — 9. are presented in Fig. [5] for 
L = 32, 64 and 128. Here the scaling collapse in the 
inset is obtained with z = 2.7, whereas the gap exponent 
is given by + 1 = 0.74, so that the localization condition 
in Eq.|j2U is satisfied. Also the Frechet distribution fits 
very well the scaled curves. 



3. Random Heisenberg models in Id and 2d 

The random quantum models studied so far have a dis- 
crete symmetry. In this section we are going to consider 
random Heisenberg models which have continuous sym- 
metry and the renormalization procedure and the cor- 
responding fixed points are somewhat different than for 
discrete symmetry models. To be specific the Heisenberg 
models we study here are defined by the Hamiltonian: 

Hh = 22 JijUtjSi ■ Sj , (22) 

in terms of spin-1/2 variables. Si at site, i, and the di- 
lution variables: ij = with probability, p, and t% — 1, 
otherwise. Here we consider two different types of mod- 
els, i) For the non-diluted models with p = 0, the random 
couplings are both ferromagnetic and antifcrromagnetic 
and their average is [J] av = 0. ii) For the diluted models 
we have < p < 1 and the random couplings are only 
antiferromagnetic, J^.j > 0. 

We start with non-diluted models in Id and 2d and 
calculate the smallest gap in the system by the numerical 
application of the RG method. As described in detail in 
Ref.[33| for these models the RG scales into a so called 
large spin fixed point |34| . having an effective moment, 
S e ff ~ L d< ^, with £ « 1/2. Due to the formation of large 
spins the low-temperature singularities of these models 
are also different than that of the RTIM, for example the 
average susceptibility has a Curie-like behavior. 

The distribution of the log-gaps for the Id non-diluted 
model is shown in Fig. |S|for different sizes, and in the 
inset the scaled curves are presented with z — 3.3. There 
are considerable deviations from the Frechet distribution, 
in particular for small gaps. The relation in Eq.J2U is 
not satisfied, so that the excitations are non localized. 

The distribution of the log-gaps for the 2d system is 
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FIG. 6: (Color online) Distribution of the log gaps, obtained 
at the last step of the RG procedure for the Id non-diluted 
Heisenberg model with uniform disorder, —0.5 < Ji < 0.5, 
and for different sizes. In the inset the optimal scaling col- 
lapse is obtained with z = 3.3. For small gaps there is a 
considerable deviation from the Frechet distribution, shown 
by a full line. The asymptotic slope of the curves is given by 
uj + 1 = 0.67, so that the relation in Eg. 12111 is not satisfied. 

shown in Fig. 0for three sizes, L = 16, 32 and 64 and in 
the inset the scaled curves are presented with z = 2. The 
Frechet distribution seems to give a correct description 
and also the the relation in Ea. i|21|) is satisfied. We note 
that in 2d there is frustration in this model, thus it is 
called as a spin glass. 
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FIG. 7: (Color online) The same as in Fig|^]for the square 
lattice Heisenberg spin glass with Gaussian disorder of vari- 
ance 1 for different sizes. In the inset the optimal scaling 
collapse is obtained with z = 2. The Frechet distribution is 
indicated by a full line and has a satisfactory fit to the data. 
The asymptotic slope of the curves is given by w + 1 = 1, so 
that the relation in Eq, 12H is satisfied. 

Next, we consider diluted models on the square lattice 
for which the couplings are antifcrromagnetic and dis- 
tributed uniformly, < < 1. For the dilution, p, we 
consider two values. In Fig|S]we present the distribution 



of the gaps at p = 0.125, which is below the percola- 
tion threshold, whereas in FigElwe consider the dilution 
at p = 0.45, above the percolation threshold, when the 
system is broken into non-interacting finite clusters. 

In the first case during renormalization the system 
scales into a large spin fixed point. As seen in FiglHlthe 
distribution of the calculated gaps differs considerably 
from the Frechet distribution and also the the relation in 
Eq.(j2U is not correctly satisfied. 
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FIG. 8: (Color online) The same as in FigE] for the square 
lattice Heisenberg antiferromagnet with dilution p = 0.125. 
In the inset the scaling collapse indicates z = 5. The Frechet 
distribution, shown by a full line, differs considerably from the 
numerical data. The asymptotic slope of the curves is given 
by u + 1 = 0.35, so that the relation in Eq. l21H is satisfied 
only with some error. 

For p = 0.45 the system is separated into indepen- 
dent parts and the smallest gap of the system is just 
the smallest gap of these clusters. Since the energy gaps 
of the clusters are expected to be distributed in identical 
power-law form, the applicability of the EVS is probable. 
Indeed, in Fig^Jthe numerically calculated gap distribu- 
tions can be well described by the Frechet form and also 
Eq.J^U is satisfied. 



III. GRIFFITHS SINGULARITIES IN RANDOM 
STOCHASTIC SYSTEMS 

A. PASEP with extreme disorder 

First, we consider the PASEP with pw disorder as de- 
fined above Eq.Q with a special form of the bimodal 
disorder, when the particles are of two kinds. For a frac- 
tion of c <C 1 black particles the hopping rates are: pi = 1 
and qi = A, whereas the (1 — c) fraction of white particles 
have hopping rates: pi = 1 and qi = A -1 . In order to ob- 
tain exact results we take the limit A> 1. With this type 
of disorder the drift of the white particles to the right is 
slowed down by the black ones and the slowing dawn is 
even more effective if two or more black particles happen 
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FIG. 9: (Color online) The same as in Fig|H| for dilution 
p = 0.45, above the percolation threshold. In the inset the 
optimal scaling collapse is obtained with z — 10. The Frechet 
distributions, indicated by full lines, fit very well the numer- 
ical data. The asymptotic slope of the curves is given by 
uj + 1 — 0.2, so that the relation in Eq. l2H is satisfied. 



to stay behind each other. A cluster of n black particles 
has a very small density: p(n) — c n and its speed to the 
right, in the absence of another black particles and clus- 
ters is given by, v(n) rj A -t *. Indeed in order to make 
the complete black cluster one step to the right all the n 
particles should perform their hop to the unpreferential 
direction at the same timestep. At this point it is easy to 
notice the isomorphism of this problem with the RTIM 
in Sec lII A H with the correspondence: e <-> v. Thus we 
can immediately write for the distribution of the effective 
velocities of the black clusters: 



P(v) = ±f 







(23) 



with the exponent, oj, defined in Eq.JSJ). 

Evidently, the stationary velocity of the PASEP is 
given by the smallest effective speed of the clusters, v\, 
and all the particles move behind that slowest black clus- 
ter. In a large finite system with a finite density of 
particles, the stationary velocity goes to zero as given 
in Eq.||BJ with a dynamical exponent defined in Eq.@. 
Consequently in a finite system the distribution of v\ is 
given in the scaling form in Eq. I|1U|) , and the distribution 
function of the scaling variable: u = uqV\L Zv is given by 
the Frechet distribution in Eq. flTJl , 

The results obtained in this section can be easily gen- 
eralized for the PASEP with site- wise (sw) disorder, in 
which case the hop rates depend on the position: for 
the site i they are pi (right) and tfo (left). The control- 
parameter of the model is the same as for pw disorder in 
Eq. |@J . For the extreme binary disorder considered above 
most of the sites are white and promote the movement 
of the particles to the right, but the few black sites and 
in particular the rare black clusters form barriers, which 
slow down the particle motion. This problem is studied 



in more details in Ref. t 9j . Here we just note that the aver- 
age velocity of a particle which goes trough a single large 
barrier of size, n, is given byj^, v w A"™/ 2 , since due to 
particle-hole symmetry in the stationary state the barrier 
is filled up to n/2, so that the particles should make n/2 
consecutive steps against the barrier. As a consequence 
the derivation in the previous paragraphs for pw disorder 
should by slightly modified, which leads to a dynamical 
exponent, z s , given by: 



~2~ 



(24) 



In particular the distribution of the stationary velocity in 
a finite system is still given by the Frechet form in terms 
of the scaling variable: u = uqV\L Z3 . 



B. Strong disorder RG and scaling results 

Here we show that the results in the previous section, 
i.e the relation between Griffiths singularities and EVS 
holds for a general form of disorder, too. We start with 
the PASEP with pw disorder for which a variant of the 
strong disorder RG approach has been applied 8]. For 
this model during renormalization the fastest hop rates 
are consecutively decimated out and new clusters of par- 
ticles are created with effective hop rates obtained by 
a perturbation calculation. Without going to the details 
we mention that there is a one-to-one correspondence be- 
tween the RG rules for the RTIM and that of the PASEP. 
In the Griffiths phase for S p > in the first part of the 
renormalization both left and right hop rates are deci- 
mated out until particle clusters of typical size, £, are 
created. (Here £ is the average correlation length, which 
behaves for small S p as £ ~ S~ 2 .) At this scale of the RG 
the remaining effective particles have a practically van- 
ishing left hop rate, qi « 0, and finite effective right hop 
rates, which have an asymptotic power-law distribution: 
P(p) ~ p- 1 + 1 / z p. Here the dynamical exponent, z p , is 
given through the equation: 



l/Zp 



(25) 



which makes the analogy with the RTIM complete, see 
in Ea. H14fl . The renormalized PASEP than consists of in- 
dependent particles the velocity (right hop rate) of them 
has the same power-law distribution and the smallest of 
them is the stationary velocity of the finite system. Con- 
sequently the conditions of EVS are asymptotically sat- 
isfied, so that the distribution of the stationary velocity 
in a finite system is given by the Frechet form in terms 
of the scaling variable: u — Mo^i(i/£) 2p . 

For the PASEP with sw disorder one can not sim- 
ply apply the strong disorder RG, so that we use here 
phenomenological, scaling considerations. As noticed in 
Sec IHI Al the rare regions are represented by large bar- 
riers, which in the length-scale, £, are expected to be 
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independent and well separated from each other. Scal- 
ing consideration in Ref.0 show that the distribution 
of the velocities associated to large barriers is given 
by, P(v) ~ ^-i+i/zs anc j f or dynamical exponent, 
z s , the relation in Ea. l24() holds in the entire Griffiths 
phase. Once more the stationary velocity in a finite 
system is given by the smallest velocity associated to 
the largest barrier, and its distribution is expected to 
be in the Frechet form in terms of the scaling variable: 
u = UoVtiL/^Y". 

Finally, we consider the Id zero-range process (ZRP) 
with quenched disorder 10J . In this model the i-th lat- 
tice site can be occupied by rii > particles, from which 
the topmost one can hop to nearest neighbor sites with a 
position dependent rate: qi to site i + 1 and to site 
i — 1. It is known that the ZRP with this definition can 
be exactly mapped (up to translations of the configura- 
tions of the lattice) to the PASEP with pw disorder, as 
studied here. The sites of the ZRP are particles in the 
ASEP and the particle clusters in the ZRP correspond 
to holes in front of the particles in the PASEP. Then the 
stationary current in the ZRP is just the stationary par- 
ticle velocity of the PASEP. From this mapping and the 
previous reasoning follows that the stationary current of 
the ZRP, J, in a finite system scales as: J ~ L~ z p, and 
its distribution is the Frechet distribution. 



C. Numerical results 

We start to analyze the distribution of the station- 
ary velocity, v, of the PASEP with uniform pw disor- 
der, which is presented in Fig^| Here we have evalu- 
ated an exact expression for v, which is described c.f. in 
Ref. @ and in this way we have studied chains with up to 
N = 2048 particles over 100000 independent realizations 
of the disorder. Fig^|shows an excellent agreement with 
the Frechet distribution having the exact dynamical ex- 
ponent in Eq. (|23l) . 

For the PASEP with sw disorder there is no analytical 
expression for the stationary velocity so that v is calcu- 
lated by simulation. We have considered 10000 random 
half-filled chains with binary disorder of different lengths 
up to L — 512. The results as presented in Fig^Jare in 
good agreement with the Frechet distribution, in which 
the exact dynamical exponent is taken from the scaling 
result in Eqs.|22|) and CSJ- 



IV. DISCUSSION 

In this paper we have considered strong Griffiths ef- 
fects in different interacting many particle systems and 
studied their possible relation with extreme value statis- 
tics. Our examples included random quantum systems 
in one and two dimensions as well as stochastic systems 
with quenched disorder in one dimension. Our exact, nu- 
merical and RG results indicate that for systems having 




-4 -2 2 4 6 8 
x=-ln[vN z ] 



FIG. 10: Distribution of the stationary velocity of the PASEP 
with pw uniform disorder having ho = 3. calculated for differ- 
ent numbers of particles. The Frechet distribution with the 
exact dynamical exponent in Eg. 12511 is indicated by the full 
line. 
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x=-ln[vL z ] 



FIG. 11: The same as in FigEE for the PASEP sw binary 
disorder with c = 0.25 and A = 2. For the Frechet distribution, 
which is indicated by the full line the dynamical exponent is 
taken from the scaling result in Eqs. 1241 and 1251 . 



a discrete symmetry the distribution of the inverse time 
scales (excitation energy for quantum systems, station- 
ary velocity for exclusion models) in large finite samples 
has a universal form, which is the limit distribution of the 
extremes of iid random numbers. In these examples the 
Griffiths singularities are characterized by the dynami- 
cal exponent, z, which is a continuous function of the 
control parameter, S. However the distribution function 
depends on, u — L z / d T-\ and given by the Frechet dis- 
tribution in Eq. (|llf> . The physical picture behind this 
result is given by the strong disorder RG method: dur- 
ing renormalization fast degrees of freedom are gradually 
decimated out and the system finally transforms into a 
set of practically independent degrees of freedom. The 
characteristic time-scales of these localized units follow 
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a power-law distribution with a z-dependent decay ex- 
ponent. Since the physically relevant relaxation time is 
given by the largest time-scale we arrive to the results of 
EVS. 

The universality of the distribution function of differ- 
ent problems is tested by numerical and RG calculations. 
For models with a discrete symmetry in all cases we ob- 
tained convincing evidence of universality. On the con- 
trary for the random Heisenberg model, which has a con- 
tinuous symmetry, the distribution function is found to 
depend on the specific form of the disorder. Interestingly, 
for the 2d spin-glass model as well as for the diluted 2d 
model above the percolation threshold the distribution 
function is found in universal Frechet form. 

One might ask the question how general these results 
are. On the basis of the RG approach we conjecture that 
for all interacting systems which have a disordered Grif- 
fiths phase the singular properties of which are controlled 
by the same type of strong disorder fixed points as for 
the RTIM the distribution function of the inverse time- 
scales is universal. Possible systems of this class are, be- 
sides random quantum magnets and exclusion processes, 
some reaction-diffusion models the dynamics of the 
random-field Ising chain 36] , the localization of a random 
polymer at an interface[23, etc. 

Evidently the above considerations of EVS does not 
apply for the distribution of average and local physical 



quantities in random systems. For the RTIM average 
quantities are, am ong others the uniform susceptibility or 
the sound velocity whereas examples for local quan- 
tities are the surface susceptibility or the surface magne- 
tization in the ordered Griffiths phase |26|, l39l l40l . The 
biased random walk in a random environment [4 lj is a 
one particle problem, thus the applicability of the EVS 
is not expected. Indeed, the stationary velocity is an av- 
erage quantity, since the time needed for the particle to 
get through a system is obtained by averaging the wait- 
ing times associated at different points of the lattice[42|. 
Therefore the distribution of the stationary velocity is 
not in the Frechet form[4l|. 
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